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Abstract 

Nonlinear diffusion dtp = A(<f>(/3)) is considered for a class of nonlinearities 
It is shown that for suitable choices of an associated Lyapunov functional can 
be interpreted as thermodynamic entropy. This information is used to derive an 
associated metric, here called thermodynamic metric. The analysis is confined to 
nonlinear diffusion obtainable as hydrodynamic limit of a zero range process. The 
thermodynamic setting is linked to a large deviation principle for the underlying 
zero range process and the corresponding equation of fluctuating hydrodynamics. 
For the latter connections, the thermodynamic metric plays a central role. 


1 Introduction 

There is a wealth of results on nonlinear diffusion equations of the kind 

Ap = A(4>(p)); (1) 

as the monumental monograph by Vazquez m for the porous medium equation, 
$(p) = p™, demonstrates. Recently, Bodineau, Lebowitz, Mouhot and Villani |1] 
studied drift-driven nonlinear diffusion equations with inhomogeneous Dirichlet data 
and established the existence of Lyapunov functionals for such equations. 
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Starting with the work of Jordan, Kinderlehrer and Otto [20], diffnsion eqnations 
have been written as gradient flow of the entropy. For the example of linear diffnsion, 
dtp = Ap, with the Boltzmann entropy density Fun := plog(p), it is immediate that 
dtp = div(pVS'(;^). The seminal insight of [20] is that the operator involved, K{p)^ := 
—div(pV.^), is associated to a metric and hence dehnes a geometry; the evolution is 
then determined by the steepest descent of the entropy S in this associated geometry. 
Here S is the classic entropy given by iS(p) = / S\i-a_{p{x))dx if p is absolutely continuous 
with respect to the Lebesgue measure and S{p) = +oo if p is not absolutely continuous. 
More generally, gradient flows are formally equations of the form dtp = —K{p)dS{p), 
where iS is a functional dehned on the state space, dS the Frechet differential and K (p) 
an operator from the cotangent space to the tangent, the so-called Onsager operator. 
This formulation says more than that iS is a Lyapunov functional; it says the evolution 
is given by the steepest descent evolution of S in the geometry given by K. 

In this short note, we ask a simple question: For equation ([1]), with a suitable 
choice of $, can one show that an associate Lyapunov functional is in fact a thermody¬ 
namic entropy, that is, the entropy associated with an underlying particle system? This 
requires that ([1]) is the hydrodynamic limit of the underlying particle system. Further¬ 
more, we will then write ([T]) as formal gradient flow (steepest descent) of the entropy, 
in a metric we derive as the metric associated with the thermodynamic entropy. We 
thus call this metric (or geometry) the thermodynamic metric (or geometry). In some 
cases, the formal setting can be shown to dehne rigorously a gradient flow. 

We employ basic tools, techniques and results from statistical mechanics and op¬ 
timal transport. Yet, the result seems interesting in a twofold way. Firstly, from the 
statistical mechanics perspective, the thermodynamic entropy is well known (see, e.g., 
Grosskinsky m), but has not been placed in a gradient flow context, i.e., linked to the 
associated thermodynamic metric. Secondly, the nonlinear diffusion dtp = A$(p) has 
been studied in a classic Wasserstein setting, with the 2-Wasserstein metric [2], and 
has been shown to be the gradient flow, with respect to the 2-Wasserstein metric, of 
the free energy given by 

Hp)^ j F(p)d^ + F'(^)ppT‘). 

Here F is characterised by the ordinary differential equation $(p) = pF'{p) +F{p) and 
assumed to satisfy F'(oo) := sup^^o (^^e [2l Theorem 11.2.5] for less restric¬ 

tive assumptions). This gradient flow setting is different from the setting presented in 
this article. 

It is well known that one equation can have different incarnations as gradient flow, 
that is, as steepest descent of a functional in a metric space; linear diffusion is a classic 
example with inhnitely many such formulations. Then the functional is a Lyapunov 
functional. Here we show that in the setting of this article, the Lyapunov function has 
a thermodynamic interpretation, being the thermodynamic entropy of the so-called 
zero range process. We point out that the existence theory can be established without 
this thermodynamic interpretation. That is, the knowledge of Lyapunov functionals. 
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which may or may not be entropies, can be key to a successful existence theory, as [13] 
testihes. At the same time, large deviation theory singles out one particular metric in 
the cases under consideration. We show in Section 0] that this gives, at least formally, 
immediately rise to a variational scheme of Jordan-Kinderlehrer-Otto (JKO) type [2Uj . 

We remark that different particle processes can lead to the same equation as hy¬ 
drodynamic limit, but that equation will then typically have different representations 
as pairs of metric and entropy. For example, the symmetric exclusion process and a 
system of independent random walks both lead to the linear diffusion equation, but 
with two different entropies and correspondingly two different metrics, see [T]. 

It turns out that in the setting of this article, the same pair of a weighted Wasserstein 
metric and the thermodynamic entropy arises both from large deviation theory and a 
weighted version of the classical Benamou-Brenier formula. The weighted Wasserstein 
metric is (formally) conformally equivalent to the classical one, with the weight being 
related to the diffusion of a tagged particle. Namely, it is often convenient to write ([I]) 
as 

dtp = A(a(p)p) = div(a(p)Vp) + div(pV ct(p)), (2) 

since a{p) = <h(p)/p will be related to the diffusivity of a tagged particle in the underly¬ 
ing zero range process, as discussed in Subsection 12.21 Roughly speaking, the equation 
which is the hydrodynamic limit of the zero range process looks like the Kolmogorov 
equation of a single tagged particle. This diffusivity (sometimes called self-diffusion co¬ 
efficient) will turn out to be closely related to the thermodynamic metric of the process 
discussed below. 

Plan of the paper We focus on the so-called zero-range process as underlying micro¬ 
scopic model; this process and its properties are described in Section [2l In particular, 
the entropy (density) S is identihed there; in Section [3] the associated metric is deter¬ 
mined. We show that this metric can also be obtained differently, by a large deviation 
argument; this is the content of Section 01 which also sketches connections to associated 
models of fluctuating hydrodynamics. 

2 Particle models 

2.1 The zero range process 

In this subsection, we briefly summarise the zero range process on a flat torus; see m 
[TB] for more information on zero range process and the condensation phenomenon they 
exhibit for particular choices of jump rates. 

A zero range process is an interacting particle system. Let A be lattice, that is, an 
ordered collection of boxes indexed by an (at most countable) index set. Each box can 
be empty or contain a hnite number of particles. The particles are indistinguishable. 
They evolve according to an irreducible Markov kernel on A in the following way. A 
particle at a given site x G A interacts only with the particles located at the same site 
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(hence zero range). Specifically, the rate at which one particle leaves a site x depends 
only on the total number of particles at site x, described by the (local) jump rate 
function 5^: Nq := N U {0} — > Rjj', which satisfies g{k) = 0 if and only if /c = 0. Here 
g is assumed to be Lipschitz continuous, i.e., sup^^pj^^ \g{k + 1) — g{k)\ < +cxo. 

In particular, for the hydrodynamic limit in the fiat torus we consider the diffu¬ 
sively rescaled zero range processes with finite lattice the discrete d-dimensional torus 
^x)'^ = {0,1,..., iV — 1}'^ as N ^ - 1 - 00 . The state space of the system is 

T'* j 

the set Mtv := Nq^ of all configurations of particles 77 : —> No; the system is in 

state ?7 if x contains ri{x) particles for all x G T^. Let p denote the symmetric nearest 
neighbour transition probability (renormalised so that it has total mass equal to 2d) 


d 

pix,y) :='^l{-ej,ej}{y-x). 
i=i 


The symmetric (nearest neighbour) zero range process with jump rate g is the Markov 
jump process on the state space with formal generator 

LNfiv) = - f{v)}9{v{x))p{y - x), 


where /: M^v —> M is a function on the space of configurations, 77 ^’^ := 77 — Ij^.} + l^yj 
is the configuration resulting from 77 by moving a particle from x to y and the term 
comes from the diffusive rescaling of the process. 

The microscopic stochastic dynamics defined by this generator can be described as 
follows. Starting from an initial configuration 770 , the system waits an exponential time 
Ti with rate := which time one particle is moved from 

X to y with probability ^^j^^yp{y — x). Then the processes starts again with 77.^1 := 
in place of 770 . 

The equilibrium states of zero range processes are well known. They are charac¬ 
terised by the equation uLn = 0 Prop. A.4.1]. If 1 / is a translation invariant 
measure on the configuration space Mat, all its one-site marginals r'xi') ■= ^{v{x) = •}, 
X G T^, must be equal. In particular, ip := W,^g{ri{0)) = K^g{ri{x)). So if z/ is a trans¬ 
lation invariant and product equilibrium state with fugacity (p := W.^g{ri{0)) , then 
uLn = 0 can be solved to yield 


z/{77 77(0) 


k} = 7/(77 


77(0) = 0 } 


d!(fc)’ 


(x, k) X No 


where g\{k) := g{l) • ... 
partition function 


■ g{k) and ( 71 ( 0 ) = 1. Since 

00 

z('f) ■= Y. 


00 1. 


k=0 




z/ is a probability measure, the 
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converges at and ^{ri I ri{0) = 0} = 


z{v)- 


Conversely, for ip in the proper domain 


Dz ■= {</? G Kq I Z{cp) < +oo}, the translation invariant product measure 
with one-site marginal 


N ■ 


m 


\1N 


{^( 0 ) = k} := 




Z{ip)g\{k) 


is an equilibrium state of the zero range process. Consequently the family of the trans¬ 
lation invariant and product equilibrium states of the zero range process is {k^}^£Dz- 
From a statistical mechanics point of view, this family is the grand canonical ensemble. 
Of course, in order for non-trivial equilibrium states to exist, the critical fugacity 


(Pc := liminf \/g\{k) 

k^+oo 


(3) 


must be strictly positive to ensure that the partition function Z has non-trivial radius 
of convergence. 

For the description of the evolution of the empirical density of the zero range pro¬ 
cess, it is useful to have a reparametrisation of the grand canonical ensemble by the 
density. This can be achieved by introducing the density function R : Dz —> [0, -|-oo], 
which is given by 


:= Ep^r7(0) 


Z{p>) 


oo 




gm 


g^z'{ip) 
Z((p) • 


(4) 


This function is analytic (on [0, p>f)) and strictly increasing [211 Section 2.3] with Dr = 
Dz n Dz'- The critical density is then := sup R{Dr) G (0,-|-oo]. In this article, 
we restrict our attention to zero range processes with superlinear jump rate g, i.e., for 
some constant oq > 0 

g{k) > aok for all k > 1, (5) 

which enforces p>c = Pc = +oo (we refer the reader to [16] for zero range processes that 
exhibit phase transitions). Then the inverse <h := R~^ is well dehned on all of M||; by 
setting Pp := the grand canonical ensemble is reparametrised by the density. We 
note that $(p) is stochastically represented as the mean local jump rate with respect 
to the translation invariant and product equilibrium state of density p, since 






z(4W) tr gm 


( 6 ) 


The thermodynamic entropy density S of the zero range process is dehned as the 
Legendre transform of the pressure P{\) = log(Z(exp(A)); here it is (see [161 Eq. (3.2)]; 
we use the opposite sign convention) 


S{p) = plog$(p) - log(Z(<h(p))) (7) 


5 










and thus 


S'(p)=log(4W)+p||d_|(|(|4-(,). 

Since $ = -R“\ inserting ip = <l>(p) in (j4]) we get 


( 8 ) 


and so 


and hence (0) simplihes to 




Z(4(p)) 


4>(p)’ 


S'(p) = log(<l>(p)) = log(cr(p)p). 


(9) 


The large deviations functional with speed w (see also Subsection 14. 1 1I of the 
occupation variables {ri{x )}with respect to the equilibrium state on the inhnite 
lattice is given by [2T1 Lemma 6.1.7] 


IpSp) = Plog 


^(P) 


- log 


zmp)) 

z{^{p*)r 


Thus it differs from the thermodynamic entropy by an affine function that depends 
only on the choice of the equilibrium state through its density p*. 

The integral functional S corresponding to the thermodynamic entropy S is given 
by 

S{p)= ! S{pndx + S\^)p^{T^), (10) 

Jld 

where S\oo) = hmp_j.+oo S{p)/p = log<pc is the rate of linear growth at infinity of the 
thermodynamic entropy, and p = p“ + p^ is the Radon-Nikodym decomposition of 
p, p“'^ <C Cfd, p^ T Cfd (with respect to Lebesgue measure on the torus). In 
particular, in the case that <pc = +oo, the thermodynamic entropy has superlinear 
growth S\oo) = +CX), and thus S{p) = +cx) if p ^ Cjd. The integral functional Xp^ 
associated to the rate functional Jp, as in (ITOD . 


Xp,{p)= ! /p.(p“)da; + /;((^)p^(T'^), 

is (see the proof of Lemma 5.1.6 in [21]) the rate functional for the large deviations 
principle (with speed ^) satisfied by the sequence 

(4 E eP(X+(T'')). Neti, 
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of the empirical embeddings of the zero range process with respect to an eqnilibrium 
state of density p*. In the last display Ai + (T‘^) denotes the set of all non-negative 
Borel measures on the torus, and V the set of probability measures; the push-forward 
/jjP of a measure p by a (measurable) map / is dehned as the measure satisfying 
= p(/“^(A)) for all (measurable) sets A. Since S and Ip^ differ by an affine 
function, the integral functionals S and Xp^ differ by a constant and consequently S 
and Xp^, p* G Mg dehne the same gradient flow (with respect to any metric). 

It is well known (see [211 Theorem 5.1.1] for a precise formulation) that the hydro- 
dynamic limit of the zero range process is given by ([T]), with <h being the mean jump 
rate function. This means that starting the process from a sequence {p^} of initial 
distributions associated to a prohle po G L^(T'^), i.e., from a sequence of initial laws 
, N >1, such that for all 5 > 0 and G G C(T'^) 


lim p^ 

V^+cx) ^ 


1 






ri{x) 


Ijd 


G{x)po{x) dx 


> (5 


0 , 


then the law pf^ of the process at each time f > 0 is associated to the prohle pt, where 
{pt)t£R+ is the unique weak solution of ([T]) starting from pg. Furthermore, if the initial 
prohle is of class for some 6* > 0, then the relative entropy method shows that 
asymptotically as N ^ -|-cxo the law pf^ looks like the product measure ^pp.) with 
slowly varying prohle pt, which is characterised by the marginals 




1 <^{pt{x/N))’^ 

Z^ptix/N))) ^!(fc) ’ 


see m Theorem 6.1.1] for a precise formulation. 


2.2 Asymptotic behaviour of tagged particles 


Given a symmetric zero-range process, we can “tag” a particle that is originally at the 
origin, say, and follow its movements, denoting its position at time t by X{t). It is 
natural to ask whether there is a central limit theorem for this tagged particle, that 
is, whether the law of N~^X{N^t) converges to a dihusion dX = a{X)dW, where W 
is a standard Brownian motion, and if so, what the dihusivity a is. This central limit 
theorem has been shown to hold in many situations, both in equilibrium (i.e., with the 
initial distribution of the particles given by an invariant measure for the process) and 
in non-equilibrium, see, e.g., [21] and the references therein. It turns out [19] that for 

p > 0, 


(^ip) 


Hp) 

P 


where <F and p{t,x) as in ([T]). In the next section, this intrinsic diffusivity will be 
related to what we will call below the thermodynamic metric of the process. 
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3 Thermodynamic formulation 

We now phrase (Ej) with nonlinear a in a Wasserstein setting. We extend the classic 
Benamon-Brenier formulation j3] of optimal transport by a weighted version. 

As a starting point, note that with the notation from the previous section (see, 
e.g., &) 

dtp = A($(p)) = div(<h(p)V^'(p)) = div(paV^'(p)) (11) 

with cr(p) = ^{p)/p. In order to simplify notation, we will often omit the dependence 
of a on p. 

In this section we will show that this structure allows us to interpret (Ej) as a 
gradient flow of the entropy flTOj) with respect to a metric which is associated to the 
operator div($(p)V-) appearing in the equation above. This metric will be shown to 
be a weighted Wasserstein metric with weight a~^. 

Here we will assume that our spatial domain is the d-dimensional flat torus For 
general bounded domains and homogeneous Dirichlet boundary data the reasoning is 
analogous. The case of inhomogeneous boundary data is covered in [1]. 

We introduce a metric tensor, following the work of Otto [23]. See also the presenta¬ 
tion in [25l Section 8]. Here we work with the set ParT'^ of probability measures on the 
torus with continuous Lebesgue densities. What follows can be seen as the Benamon- 
Brenier formulation of optimal transport with a nonconstant metric. In this context, a 
point p G ParT^ in this set has formally as tangent space the space of all infinitesimal 
variations at f = 0 of smooth curves (—e, e) 3 t pt with po = p. By the conserva¬ 
tion of mass, these variations satisfy j dt\t=opt^x = 0 and thus the tangent space as 
a set can be identified with the set of all continuous functions C ,: —)■ M such that 
f ( dx = 0. The tangent spaces are thus all equal as sets, while the metric we are going 
to define depends on the base point p. As in the standard case (e.g., [23]), a metric is 
dehned on each tangent space by associating to each tangent vector ( = dtjt=oPt (i-e., a 
derivative of a path in the space of densities/measures) a vector held representing this 
variation, via the equation 

c + div(pn) = 0, (12) 

with unknown the vector held n: ^ To ensure uniqueness of the vector held v 

representing the tangent vector (, we minimise the weighted kinetic energy under the 
constraint oa. 

inf I / pa~^ \v\^ da; I C + div(pn) = 0 I. (13) 

^ Ut<^ J 

In the classic case (no division by a), this term has an interpretation as twice the 
kinetic energy, and we thus call ffT3|) (twice) the weighted kinetic energy. We take outer 
variations of the form v -|- ep~^w with divergence-free w] note that the latter property 
ensures that flT^ is preserved. Then the extremality condition for flT3|) reads 


for all divergence-free w, where (•, •) is the usual Euclidean scalar product on W^. By 
the Helmholtz-Hodge decomposition of vector helds into divergence-free and gradient 
helds, this implies 

a-h = (14) 

for some real-valued function or equivalently v = crV^. Note that thus ( = 
—div(pt>) = —div((jpV.^). As (for given suitable () the equation —div(pcrV.^) = ( 
has a unique solution, there is (up to additive constants) only one such ^ where 
the latter notation indicates the dependence on (. In this way, we can associate to 
each tangent vector ( the vector held := cr(p)V.^(. Then on the tangent space at p, 
we dehne the metric 

9p{CiX2) ■■= [ -^{v^„vi;^)dx= f pcT(p)(V^Ci> (15) 

Jjd cyp) Jjd 

For each p, one can think of a~^{p){v,v) as a Riemannian metric on a manifold which 
is topologically This penalises transport in regions with slow diffusion, where a{p) 
is small. In the classic setting [23], the division by a is absent. As in classic formal 
Otto calculus, the cotangent space is a space of functions modulo constants and thus 
the function corresponding to a tangent vector ( via (1141) is the cotangent vector 
corresponding to ( via the metric g. 

Let now S be any functional of the form flTOl) for some integrand S : —)■ R, 

which, by the regularity assumptions made earlier on the probability measures under 
consideration, is in fact given by 

S{p)-.= f S{p)dx. (16) 

Jn 

We recall that in the usual pairing, the differential dS is represented by the vari¬ 
ational derivative = S'{p). Analogously, if we denote by the dual pairing 
between tangent and cotangent space at p, and if ( is represented by V(^ = crV^i^, an 
integration by parts yields 

{dS{p)\C)p = J C^^^dx = - j diY{pv^)S'{p) da; 

= J p{vo^S'{p))dx = J -^{vi;,a{p)VS'{p))dx 

= 9p{C, -dw{paVS'{p))). (17) 

The last equality follows by dehnition of the metric g, since if ( := —div^paVS'^p)) 
then = (t(p)VS"(p). In words, formula (1T7|) says that the cotangent vector dS{p) is 
represented via the metric tensor g by the tangent vector —div(crpVS"(p)). Thus in 
our case the Onsager operator K = corresponding to the weighted metric g maps 
dS to the tangent vector ( = —div(crpV[S''(p)]), that is 

K{p)dS{p) = -diY{apVS'{p)). 


9 






Note that the identihcation C ^ given via the minimisation of ffl^ and ffT^ is in 
fact the inverse of the Onsager operator K. 

A gradient flow with respect to a metric is distinguished by the fact that the velocity 
dipt of a path t ^ pt is related via the metric to the rate of change of a functional, 
i.e., here ( = —dw{ap\/S'{p)), and so the gradient flow equation for a functional of 
the form ffTOj) becomes 


dtpt = -K{p)dS{p) = div((TpVS"(p)). 


So now if we choose the thermodynamic entropy ([7]) as the integrand in the functional 
S, with its derivative given by ([8]), then we recover ([2]). In summary if we choose 
the thermodynamic entropy, then precisely the choice (fT^ for the metric recovers the 
nonlinear diffusion equation ([1]); we thus call (1151) the thermodynamic metric. 

Weighted kinetic energies as (IT^ have been introduced and studied before, in par¬ 
ticular by Dolbeault, Nazaret and Savare m, and Carrillo, Lisini and Slepcev [5]. We 
sketch how their results can be used to obtain existence and uniqueness for solutions 
to the gradient flow corresponding to ([T]) in the thermodynamic setting. The global 
metric corresponding to the kinetic energy (fl^ is given by 

lT^. 2 (/i, zz) := inf f f pa~^\vf dxdt, 

Jo Jt‘^ 

where the inhmum is taken over all pairs {pt,Vt)o<t<i with po = p, pi = u that satisfy 
the usual continuity equation dtp + div(pn) = 0. Equivalently, one can consider 


IT'^. 2 (p., z/) := inf / / <I>(p) |n|^ dx df 

Jo Jj‘i 

with the constraint dtp + div($(p)n) = 0; the latter variational formulation is [IHl Eq. 
(1.6)] in our notation. Then, as proved in [TOl [5], concavity of d) and a type of gener¬ 
alised McCann conditions are sufficient conditions to establish existence and uniqueness 
of the gradient flow of the thermodynamic entropy with respect to the weighted Wasser- 
stein metric hW. These generalised d-dimensional McCann conditions 0 Dehnition 4.5] 
take in our context the form 

<h(p)<h'(p) > (1 — 1/d) [ <l>'(r)^dr. 

Jo 

We note that in dimension d = 1 the McCann condition is equivalent to ‘h'(p) > 0 for 
all p > 0, which is true for all mean jump rate functions $ of zero range processes. 

Let us also note that not all mean jump rate functions of zero range processes are 
concave. Two simple examples where it is concave are $(p) = p and $(p) = p/(p + 1), 
which correspond to g{k) = k and g{k) = 1, respectively. Although the latter jump 
rate does not satisfy the superlinearity assumption ([5]), one can consider the jump rate 
geik) = g{k) + ek for some e > 0, which according to (E]) gives the mean jump rate 
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Figure 1: Graph of the density function ^ in the Evans model for b = 1, b = 2.5 

and b = 3.5 


$e(p) = $(p) +ep. More generally, as one can see by plots the mean jump rate function 
corresponding to the Evans model 

9b{k) = 

is concave for all parameters 5 > 0. The density function in the Evans model 

is given by [16], Eq. (4.5)] 


Rb{(p) 


(^2-^1 (2, 2, 2 + 5, ip) ^ 

(1 + b)2Fi(l, 1, 1 + fe, Lf) 


where 2-^1 (o, b, c; z) := is the hypergeometric function and 


(a)k 


a(a + 1) • ... ■ (a + /c — 1) if /c > 1 
1 if = 0 ’ 


is the rising Pochhammer symbol. Fig. [T] gives plots for Rh for some values of b. 

An example of a zero-range process with non-convex density fnnction is given by 
the Landim jump rate fnnction given by 

k if /c = 0,1 

if A: >2 



In this case, Pb^k) = k^ and the partition fnnction is Zb{(p) = 1 + X]^i ^ —■ i- + Li6(95)) 
where Li;, is the polylogarithmic fnnction. This resnlts in the density function 


Rb{(p) 


Ub-i{ip) 

1 -l- Lib(<y9) 


Fig. |2] indicates that Rb is convex for small & > 0, concave for large 6 > 0 and neither 
of the two for intermediate values. We close this section with the rather obvious remark 
that equations of porous medium type, 4’(p) = p"^, are not covered by the setting 
described here, as for the zero-range process <F'(0) is strictly positive. 
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Figure 2: Graph of the density function i?;, = ^ in the Landim model for b = 0.5, 

5 = 3 and 5 = 5 


4 Connection to large deviation theory and fluctu¬ 
ations 

4.1 Large deviation theory 

In Section [3l we used the thermodynamic entropy S from Section [2] and used it to 
identify the associated (weighted Wasserstein) metric ffTSj) : indeed, the calculation in 
section [3] shows that this metric is associated to the operator div($(p)V-) in flTT]) which 
gives the gradient flow formulation of the thermodynamic entropy. We now describe 
an alternative link between the zero-range process and this metric. The tool is Large 
Deviation Theory, which we sketch hrst for a toy model. This theory is a step beyond a 
central limit hydrodynamic limit theorem; on the most superhcial level, large deviation 
theory describes the most likely event of improbable events in terms of an exponential 
(un-)hkelihood. As explained with the toy model, the probability of stochastic events 
of is determined to be of the form P({a;'^(t) ~ exp (—here 

~ means asymptotic equivalence (precise formulations are given in any textbook on 
Large deviation Theory, for example i)- The functional I is called the rate functional 
and the central hnding of this section is that the rate functional for the zero-range 
process has as minimiser the nonlinear diffusion equation ([1]), which appears in the 
norm defined by the weighted Wasserstein metric flT^ . 

We first explain this link using a toy model of a stochastically perturbed one¬ 
dimensional gradient flow. 


dX^ = -V\X^)dt + ^/edly, 

where dW denotes a Wiener process. Heuristically, one expects that for small noise, 
e small, trajectories will tend to be close to the deterministic trajectory solving x = 
—V'(x). Large deviation theory describes the probability of observing a realisation, 
here a trajectory, by means of a rate functional /; the probability of hnding a given 
trajectory is essentially exp(—1/e I) as sketched now. For the toy model, the large 
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deviations rate functional depends on the final time T and is denoted It] it is 


/r(i,x)= / -{x + V\x)fds 

Jo ^ 

= f + ^ [ {V'{x)fds + V{x{T))-V{x{0)) 

Jo J J Jo 

with X denoting the time derivative of x. In essence, this means that 

inf It{x,x)^ , 

where the inhmum is taken over paths x = x{s) with a:(0) = x and x{T) = y. In order 
to hnd the most likely position at t = T, we have to minimise hrst over all paths such 
that a;(0) = x and x{T) = y and then over all terminal values x{T). 

Note that the hrst integral in flT^ is bounded from below by the distance d{x{0),x{T)) 
between the initial and the hnal point. Choosing T = h <C 1 and re-scaling, we see 
that we have to minimise 

^d\x{0),x{h)) + Id {x{h)) + ^ ^ {V'{x{sh))f ds. 

If d{x{0),x{h)) 3> h, then the last term is clearly of lower order, while if d{x{0),x{h)) = 
0{h) and V is smooth enough, then the last integral equals 

I j\v\x( 8 h)f ds = |v"(i(0)) + 0(h\ 

SO again the term relevant for the minimisation is of lower order. Formally this means 
that we can hnd the most likely position x{h) by minimising the functional [20] 

^d^{x{d),x{h)) + V (x(h)). (20) 

The analogous computation starting with the large deviation rate functional for the 
empirical measure of n Brownian particles (with n playing the role of the large parame¬ 
ter 1/e in the toy model; here the rate functional is dehned on function spaces) [7], then 
one obtains an equation with the same structure as (l20l) . with d being the Wasserstein 
norm, and the entropy S in place of V (see [m for the proof). Thus, one obtains 
the celebrated Jordan-Kinderlehrer-Otto functional 1201 as expansion of the associated 
large deviation rate functional. 

We now apply this reasoning in the context of nonlinear dihusion and more specif¬ 
ically the zero range process. The large deviation rate function for the zero range 
process is of the form 

[ \\dtp{s) - Ad>{p{s))\\%-i ds (21) 

Jo 


P({a;'^(0) ~ X, x''{T) ~ y}) ~ exp ( — 


( 18 ) 

(19) 
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(see [22] and references therein); here the weighted negative Sobolev norm is dehned 
as follows. Let for some non-negative bounded function w 

{u,v)rri := / w'Vu-'Vvdx, 

“ JTd. 

and denote by A~^u the solution of 

div(tcV(A“^M)) = M. (22) 


Then 


v)hi = - / uA^ vdx, 


ijd 


where the last identity used integration by parts and ([22]). 

As seen from the reasoning above, by expanding the square in 0211) . we obtain a 
metric term 

r 2 


$ 






dx 


ga{dtpt,dtpt), 


where g is as in flT^ . This follows since if V(; = crV^,^ is the vector held representing a 
tangent vector C at p via (ITTl) . then ( = —div(pn^) = —div(pcrV^f) = —div(<h(p) 
and thus A~J^)(C) = 

The null Lagrangian (the mixed term) equals 


lo X (P)) da; df, 

so it remains to compute A^^^^^(A<h(p)). As log<h(p) solves 


div($(p)V log <h(p)) = A<h(p), 


we have A^^^^^(A$(p)) = log(<h(p)), so if S'{p) = log(<h(p)) as in ([8|), we obtain 



[ dtpA^^{A^{p))dxdt 


[ S{p{T))dx- [ ^(p(0))da;, 


(23) 


and recover the thermodynamic entropy S of flTB]) . 

The remaining quadratic term is for small time intervals formally of lower order. 
Thus in summary, an expansion of the large deviation rate functional (|2T|) for nonlinear 
diffusion associated with the zero range process yields a variational formulation of 
Jordan-Otto-Kinderlehrer type, with the metric being the weighted Wasserstein metric 
considered in this article, and the functional being the thermodynamic entropy. 
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4.2 Fluctuating hydrodynamics 


The previous Subsection 14.II showed that the thermodynamic metric is related not only 
to the limit equation, but also to the probability of deviations from it. At least at a 
formal level, this connection is even stronger, as we now explain. Formally, the metric 
is given by a symmetric linear operator, which has an inverse that can be interpreted as 
covariance matrix. If we add a suitable space-time noise to the limit partial differential 
equation, then it formally has not only the same large deviations rate functionals, but 
hnite moments of the fluctuations (that is, the difference between stochastic process 
and deterministic limit) converge to the same process as for the original particle model. 
This is now explained in more detail. Consider the equations 

dp‘ = A(*(p‘))A + div( (24) 

9,p = A(<l>(p)), (25) 

where W is space-time white noise. We point out that the stochastic PDF (SPDE) fl24l) 
is not expected to have strong (i.e., pathwise) solutions. To our knowledge, the suitable 
notion of generalised solutions is an open problem. Consider, however, the fluctuations, 
i.e., the differences between solutions of fl24D and the deterministic PDF fl2SD . re-scaled 
with the scaling of the central limit theorem. As 

:=p-p^ 

is small, we can linearise around p and obtain formally that a limit (in law) of 
should solve the linear SPDF 

= A($'(p)Z°)dt + divi^y a{p^)p^W), (26) 

i.e., it is an Ornstein-Uhlenbeck process. Here is the self-diffusion coefficient of the 
zero range process, while a is the bulk-diffusion coefficient. 

This convergence statement for the ffuctuations can be stated in the following 
weaker form. Consider for smooth 7 ; —)■ M the process 

r- := £-‘/={7,p-p-), 

then by an application of Ito’s formula (see, e.g., [HI Chapter 4.5]), and the 

(formal) limit e —)■ 0 we get for any smooth F: M —)■ M that 

F{Y\t))-F{Y%0)) - f F\Y^{s)){Y\s),d>'ip{s))A^)ds 

Jo 

-i ru'(r»(p)){ff(p(p))p(p)V 7 .V 7 )* (27) 


is a martingale. 
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Note that in the last term, the expression (cr(p(s))p(s)V7, V7) can be written as 
11711^1 (^), which implies that the large deviations should be related to the metric with 

the inverse operator, i.e., ||7||^-i(p)5 as found in the previous subsection. 

We now relate the martingale (I 27 |) to the zero range process. We consider the 
empirical measure of the zero range process, i.e., for smooth 7 

(7.7) ^=7^E''(v) 

and the fluctuations 

It can be shown [laiiiiiH] for a large class of zero range processes that for any smooth 
test function, F{Y^) converge to a solution of the martingale problem flTTI) . In other 
words, the limit of the fluctuations of the zero range process is in law equal to those 
of the SPDE fl24|) . This can be seen as justihcation of the complicated noise term 
in fl2T|) : precisely this noise term is related to fluctuations of the zero range process via 
the martingale expression fl27|) . We remark that (|2^ is of fluctuating hydrodynamics 
type. In the linear case *h(p) = p, the multiplicative conservative noise in (124|) has 
been derived by Dean [S], and using quite different lines of arguments by Sturm and 
von Renesse [ 25 ]; the link between the noise and the geometry has been explained by 
a formal large deviation argument im. 

While the discussion in this paper is limited to a sub-class of zero range processes, 
it therefore does not cover all nonlinearities *h(p) but only those that occur as hydro- 
dynamic limit of such a process. Yet, we remark that a solution of fl 2 T|) for a certain 
range of e > 0 would formally yield a stochastic process whose hydrodynamic limit is 
a gradient flow of the entropy dehned by ([ 7 ]) with respect to the weighted Wasserstein 
metric, and whose large deviation rate functional is as in Subsection 14.11 The existence 
of solutions to fl24|) . however, is to our knowledge in general an open problem. 
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